#do plots for vacination example
#pdf(file=path("Vacination_forecast_plots.pdf"), height=8, width=6)
par(mfrow=c(4,2))
		manage_plot(z.0=z$calf.ratio[6,1,,], z.treat=z$calf.ratio[6,5,,], lower=calf.ci[1], upper=calf.ci[2],type="between",label= "Proportion of population juvenile", panel="A - Within", xlim=c(0,.4))
		
		manage_plot(z.0=z$bull.ratio[6,1,,], z.treat=z$bull.ratio[6,5,,], lower=adult.male.ci[1], upper=adult.male.ci[2],type="between",label= "Proportion of population yearling or adult male", panel="B - Within", xlim=c(0,.80))
		
		manage_plot(z.0=z$cow.ratio[6,1,,], z.treat=z$cow.ratio[6,5,,], lower=adult.fem.ci[1], upper=adult.fem.ci[2],type="between",label= "Proportion of population adult female", panel="C - Within",xlim=c(0,.8))
		
		manage_plot(z.0=z$n.s.total[6,1,,], z.treat=z$n.s.total[6,5,,],upper=3000,type="above",label= "Population size", panel="D - Above",xlim=c(0,15000))
		manage_plot(z.0=z$n.s.total[6,1,,], z.treat=z$n.s.total[6,5,,],lower=3500,type="below",label= "Population size", panel="E - Below",xlim=c(0,15000))
		
		manage_plot(z.0=z$pv.cow.s[6,1,,], z.treat=z$pv.cow.s[6,5,,],lower=.40,type="below",label= "Seroprevalance", panel="F - Below",xlim=c(.20,1))
		
		lower = ecdf(z$phi.s[6,5,,])(summary(z$phi.s,median)$stat[1,1]/2)
		manage_plot(z.0=z$phi.s[6,1,,], z.treat=z$phi.s[6,5,,],lower=lower,type="below",label= "Probability of infection", panel="G - Below",xlim=c(0,.4))
		
#dev.off()